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Prospect of probing dark energy using stochastic gravitational waves 
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We study the possibility of probing dark energy behaviour using gravitational wave experiments 
like LISA and Advanced LIGO. Using two popular parameterizations for dark energy equation of 
state, we show that with current sensitivities of LISA and Advanced LIGO to detect the stochastic 
gravitational waves, it is possible to probe a large section of parameter space for the dark energy 
equation of state which is allowed by present cosmological observations. 


I. INTRODUCTION 

Notwithstanding the evidence of late time acceleration 
of the Universe from a variety of cosmological observa¬ 
tions the lack of any definitive answer so far for 

the cause of this acceleration is quite troubling. That is 
why the goal for most of the current and future obser¬ 
vational efforts is to unravel this mystery. The idea is 
to study the effect of late time acceleration on different 
cosmological observables like distances, angles, volumes, 
growth of structures and many more. Subsequently mea¬ 
suring these observables with great accuracies, one can 
constrain the different dark energy or modified grav¬ 
ity models @ (the two most widely studied approach to 
explain the late time acceleration). In this approach, un¬ 
less one accurately handles the systematics involved in 
different cosmological observations, it is very difficult to 
get any meaningful constraint on dark energy models. 
After the latest results by Planck in 2015 Q combined 
with other non-cmb measurements 0,i , we now have 
the most stringent constraint on dark energy behaviour. 
This shows that the latest cmb-|-non-cmb data is not only 
consistent with simple cosmological constant (A), it is also 
consistent with other non-A behaviours and more inter¬ 
estingly with phantom behaviour @,0. As the issue is 
not settled yet, we are looking for data from future ob¬ 
servational set ups like EUCLID Q, DES 0, SKA [I^ 
etc, which may hopefully single out a particular class of 
dark energy behaviour decisively. 

Under this circumstances, it is always useful to look 
for other new observables that can be useful probes for 
dark energy. In this work, we study the prospects of 
probing dark energy models by observing the stochas¬ 
tic gravitational wave background. Gravitational waves 
provides a unique opportunity to see our Universe in a 
way that is not possible through cosmological observa¬ 
tions [l3l| . In recent times, lot of progress has been made 
in the process of detecting gravitational waves using laser 
interferometer. The Laser Interferometer Space Antenna 
(LISA) and the Advanced Laser Interferometer Grav- 
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itational Waves Observatory (Advanced LIGO) [T^ will 
provide a major breakthrough in gravitational wave as¬ 
tronomy. The chance to detect gravitational wave includ¬ 
ing stochastic gravitational wave has increased manyfold 
using these detectors. Hence it is interesting to study 
the effect of dark energy on the stochastic gravity wave 
spectrum and subsequently the possibility to probe dark 
ener gy m odels using the data from these two detectors 
(See [ia [13 for earlier studies in this direction). 

In this work, using two popular parameterizations for 
dark energy equation of state, we study the possibility 
of probing dark energy behaviour using LISA and Ad¬ 
vanced LIGO. We show that with the current sensitivities 
of these two detectors, it will be possible to probe a large 
portion of parameter space that is currently allowed by 
the present cmb-|-non-cmb data. 

Taking into account the tensorial perturbation in a flat 
Friedman-Robertson-Walker metric given by 

= a^(T)[dr^ — {Sij hij)dx’'dx^, (1) 

where hij is the tensorial gravitational field and is trans¬ 
verse and traceless. For a fixed wave vector k and a 
fixed polarization, hij is governed by the wave equation 

[I1E1I3 

+ [e - -)/rfc = 0 (2) 

a 

where hk^r) = /rfc(T)/a(r) and prime denotes derivative 
w.r.t the conformal time r. Subsequently the spectrum of 
the relic gravitational waves h(k, r) can be calculated us¬ 
ing h{k, t) = ^^/c|/ifc(T)|, where Ipi is the Planck length. 

Once h{k,T) is known, energy density Dn(fc) of the relic 
gravitational waves is given by [T^ [T^ [13 

■ (3) 

Here ko is the wave number corresponding to the 
present Hubble radius. To solve, equation (2), one has 
to assume a initial contribution for the relic gravitational 
wave background from inflation. The current data from 
temperature anisotropy in the GMB measured by Planck 
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gives an upper limit to r which is the ratio of the am¬ 
plitudes of the tensor and scalar fluctuations produced 
during inflation. Assuming the observed temperature 
anisotropy in CMB AT/T to be ~ 0.37 x 10“^ at multi¬ 
pole I = 2 [I^ (corresponding to anisotropies on Hubble 
scale), the spectrum h{k,T) is fixed at present as 

/i(fco,ro) = 0.37 X 10"^ X r, (4) 

where tq is the conformal time at present. We assume 
r = 0.1 in our calculations which is the upper limit on r 
as observed by Planck fTst ]. 

Next we need to know the scale factors a(T) at different 
epochs in the Universe to solve equation (2). For the first 
four stages of the universe, e.g inflation (I), reheating 
(RH), radiation dominated (R) and matter dominated 
(M), we assume the scale factors as [il[i3: 

a/(r) = Zi|r - -oo < r < t/, 

aRHir) = az{T - ti < t < trh, 

a-Rir) = ae{T - Te), Trh <t <Tr, 

anir) = am{T - TjnT, tr<t <tm, ( 5 ) 

where a determines the inflationary evolution. For a = 
—2, the inflationary period is exactly de-Sitter. But as we 
know that the inflationary period is not exact de-Sitter 
but a small deviation from that, we assume a = —1.95 
for our subsequent calculations. This value is consistent 
with the previous work by Zhang et al on spec¬ 

trum of gravitational waves assuming ACDM model. P 
determines the expansion during the reheating stage and 
(3^—1 [T^. The parameters tj,trh,tr and tm deter¬ 
mine the different epochs. 

To get the scale factor during the late time accelerating 
epoch, we solve the Friedmann equation. 



for time interval tm t < tq. Hq, & 
present day Hubble parameter and density parameters for 
matter and dark energy respectively. We assume a spa¬ 
tially flat Universe which is consistent with recent Planck 
result 0. w(a) is the dark energy equation of state and 
in this work we consider two most popular parameteriza- 
tions for w(a). These are the CPL [ij parameterization 
and the GCG parameterization [T^ given respectively by 

w(a) = Wo + Wa(l - a), (7) 

^ “a,- h(l-A,)a- 3 (i+“)’ 

where wo^Wa, As and a are the parameters involved in 
these parameterizations. One can show that the GCG pa¬ 
rameterization incorporates two interesting quintessence 


behaviour e.g. the thawing (1 -|- a < 0) and tracking 
behaviour (1 -|- a > 0) [l^ 

All the constants in eq. (5) will be fixed with re¬ 
spect to junction times by taking normalization a(ro) = 
Oo = 1 and then by equating a and da/dr at the 
junction times. Moreover, to fix the junction times, 
we assume a(To)/a(TM) = 1.33, a{TM)/a(TR) = 3454, 

a{TR)/a{TRH) = lO^'^, a{TRH)/a{Ti) = 300 |I3,[T3|. 


II. SOLUTION OF THE SPECTRUM 

For scale factors of power law form e.g. a(r) oc t“, the 
solution of eq. (2) has the form 


Hkir) = VkT CiJa-l/2{kT) +C2Ji/2-a{kT) , (9) 


where Ci and C 2 are two arbitrary constants. 

So, the inflationary, reheating, radiation-dominated 
and matter-dominated stages have the solutions respec¬ 
tively as 


^i'^,fl-i-i/ 2 (a;i) + A 2 J_(,g+i/ 2 )(a;i) 

BiJp^ + i/2ix2) + B2J-(f3^ + i/2){x2) 

Mfc ("t) = \/i3 GiJi/2(x 3) + C2J_i/2(a;3) , 

DiJo/2{x4) + D2J-3/2{x4) , (10) 


where xi = k\T — Ti|, X 2 = k{T — Tp), X 3 = k{T — Te) and 
X4 = k{T - Tm). 


For Inflationary stage, we can use the adiabatic initial 
condition, limfc^oo Mfe(''') which gives Ai = 

The arbitrary con¬ 
stants Hi, B 2 , Cl, C 2 , Di and D 2 can be fixed by equat¬ 
ing ^ik{T) and ^^(t) at the junction times. Moreover /3 
can be fixed using the normalization condition given by 
eqn (4). After determining all the coefficients we get /r 
and d/i/dr at tm which are given by 


..Late 


tm 


dr 


tm 


I kzi cos(fczi) — sin(fczi) 
cos(fczi)-I-fczi sin(fc 2 i) D 2 ^, 


2 1 
TT kzi 






kziDi + {1 — k^zl)D 2 cos{kzi) 

TT V L 

{k^Zi — l)Di + kziD 2 sin(fczi)|, (11) 


where zi = {tm — Tm)- Using these initial conditions we 
will solve eq. (2) numerically during the epoch tm < t < 
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FIG. 1: CPL Parameterization: Q,g vs. wg graphs: Left and right panels are at the frequencies 2.92 x 10“^ Hz (i.e. lowest 
point of the LISA sensitivity curve) and 64.56 Hz (i.e. lowest point of the advanced LIGO sensitivity curve) respectively. For 
both the graphs continuous-red, dotted-blue, dashed-black, dot-dashed-orange and long-dashed-purple lines are for Wa values 
-1, -0.5, 0, 0.5 and 1 respectively. Brown-dotted lines in left and right panels are for the amplitude of fig at lowest point of the 
LISA and advanced LIGO sensitivity curves respectively. 
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FIG. 2: GCG Parameterization: fig vs. As graphs: Left and right panels are at the frequencies 2.92 x 10~® Hz (i.e. lowest point 
of the LISA sensitivity curve) and 64.56 Hz (i.e. lowest point of the advanced LIGO sensitivity curve) respectively. For both the 
graphs continuous-red, dotted-blue, dashed-black, dot-dashed-orange and long-dashed-purple lines are for a values -0.9, -0.5, 0, 
0.5 and 0.9 respectively. Brown-dotted lines in left and right panels are for the amplitude of fig at lowest point of the LISA 
and advanced LIGO sensitivity curves respectively. 


tq using the Hubble parameter during the dark energy 
dominated epoch given by eqn (6). 


III. RESULTS 

We fix = 0.7, = 0.3, a = -1.95 and r = 0.1 

in our calculations. These values are consistent with the 
latest Planck results 0. 

In fig. (1) we show fig as a function of wq for different 
values of Wa for CPL parameterization. It depends on 
the frequency of the observed gravitational waves. In the 
left panel, we plot it for v = 2.92 x 10“^ Hz which is the 
frequency corresponding to the lowest bound for LISA 
sensitivity curve . The horizontal line corresponds to 
the lowest bound for flaw for LISA. In the right panel, we 
do the same for Advanced LIGO with the corresponding 


frequency fixed at v = 64.56 Hz [2l|. The regions above 
the horizontal line are detectable for the present sensi¬ 
tivity for LISA and Advanced LIGO. In other words, if 
LISA and Advanced LIGO measures gravitational waves 
with these sensitivity at the corresponding frequencies, 
then the regions above the horizontal line can be probed 
and constrained by gravity wave experiments. In fig. (2), 
we redo the analysis for GCG parameterization where 
we plot fig as function of for different values of a. 
We should point out that for GCG, —Ag represents the 
present value wq of the dark energy equation of state. 


Next we show how these regions occur in the current 
allowed parameter space. In fig. (3), we show this for 
GGG and CPL parameterizations. In the left panel, we 
show the current allowed region in the Ag — a parameter 
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FIG. 3: Left Panel: The allowed parameter space for GCG parameterization for JLA(SN)+BAO+CMB-Shift (dashed) and 
JLA(SN)+BAO+GMB-Shift+/cr8 (dotted) data (See [ 2 ^ for details about the data and the corresponding analysis). The 
Horizontal lines are those probed by LISA(bigger lines) and Advanced LIGO(smaller lines). Right PanehThe allowed parameter 
space for GPL parameterization for JLA(SN)+BAO+CMB-Shift (dashed) and JLA(SN)+BAO+CMB-Shift+/cr8 (dotted) data. 
The Solid outer(inner) regions are those that can be probed by LISA(Advanced LIGO) experiments. 


space for GCG for the present observational data from 
SN+BAO+CMB-Shift+/(T8 (for details about the data 
and the subsequent analysis, please see [l^)- In the same 
plot, we show the regions from fig.(2), that can be probed 
by LISA and Advanced LIGO in future. It is clear that 
for GCG, future gravity wave experiments can only probe 
thawing (1 + a < 0) region. Most of the tracker region 
(l + a > 0) allowed by the cosmological data is away from 
the region that can be probed by gravity waves except a 
very narrow region near a ^ — 1. For GPL parameter¬ 
ization, we show the same in the right panel. Here the 
regions that can be probed by LISA or Advanced LIGO 
fall right inside the region that is currently allowed by 
the cosmological data. We should stress that the GPL 
parameterization is used by almost all the cosmological 
observations to put constraints on dark energy. Hence the 
fact that the future gravity wave experiments can probe a 
considerable region in its parameter space that is allowed 
by present cosmological data is extremely encouraging. 

Although we fix r = 0.1 in our calculations which is 
the upper bound on r by Planck, but the results do not 
differ much even if we have lower values for r e.g r = 0.05. 
Moreover we fix ,0 = —1.95 in our calculations which fixes 
the cosmological evolution during inflation. For smaller 
values of the regions that can be probed by gravity 
wave experiments become smaller and for larger values 
of /3, the regions get bigger. 

IV. CONCLUSION 

To conclude, we show the possibility of constraining 
dark energy equation of state using the future gravita¬ 


tional waves experiments. Using two popular parameter- 
izations for dark energy equation of state, we show that 
detectors like LISA and Advanced LIGO can have the 
power to probe regions in the parameter space which are 
currently allowed by different cosmological observations. 
For GGG case, mostly the thawing region can be probed 
by the experiments whereas for the GPL parameteriza¬ 
tion which is the universal parameterization for dark en¬ 
ergy equation of state used by all cosmological observa¬ 
tions, future gravity wave experiments can probe a large 
section of the parameter space that is currently allowed 
by the cosmological observations. Our investigation is 
not fully robust as we fix different cosmological parame¬ 
ters to do the analysis. But probably for the first time, 
our results show the promise of probing relevant param¬ 
eter space for dark energy behaviour using gravity wave 
experiments. This definitely motivates one to do more 
detailed analysis in this regard. 
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